CONSTRUCTING PAIRS OF DUAL BANDLIMITED 
FRAME WAVELETS IN L 2 (R n ) 



JAKOB LEMVIG 

Abstract. Given a real, expansive dilation matrix we prove that any bandlimited 
function ip 6 L 2 (M. n ), for which the dilations of its Fourier transform form a partition 
of unity, generates a wavelet frame for certain translation lattices. Moreover, there 
exists a dual wavelet frame generated by a finite linear combination of dilations of %p 
with explicitly given coefficients. The result allows a simple construction procedure 
for pairs of dual wavelet frames whose generators have compact support in the Fourier 
domain and desired time localization. The construction relies on a technical condition 
on ip) and we exhibit a general class of function satisfying this condition. 



1. Introduction 

For A G GL n (W) and y G M n , we define the dilation operator on L 2 (R n ) by D^f(x) = 

1 /2 

|det A\ f{Ax) and the translation operator by T y f(x) = f[x — y). Given a n x n real, 
expansive matrix A and a lattice of the form T = PI/ 1 for P G GL n (M), we consider 
wavelet systems of the form 

{D A jT^} jeZ/YeV , 

where the Fourier transform of ip has compact support. Our aim is, for any given 
real, expansive dilation matrix A, to construct wavelet frames with good regularity 
properties and with a dual frame generator of the form 

b 

(1) </> = Y / c j D^ 

j=a 

for some explicitly given coefficients Cj G C and a, b G Z. This will generalize and 
extend the one-dimensional results on constructions of dual wavelet frames in [15, 18] 
to higher dimensions. The extension is non-trivial since it is unclear how to determine 
the translation lattice F and how to control the support of the generators in the Fourier 
domain. This will be done by considering suitable norms in IR n and non-overlapping 
packing of ellipsoids in lattice arrangements. 

The construction of redundant wavelet representations in higher dimensions is usually 
based on extension principles [7,8, 10-14, 16, 17]. By making use of extension principles 
one is restricted to considering expansive dilations A with integer coefficients. Our 
constructions work for any real, expansive dilation. Moreover, in the extension principle 
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the number of generators often increases with the smoothness of the generators. We 
will construct pairs of dual wavelet frames generated by one smooth function with good 
time localization. 

It is a well-known fact that a wavelet frame need not have dual frames with wavelet 
structure. In [20] frame wavelets with compact support and explicit analytic form are 
constructed for real dilation matrices. However, no dual frames are presented for these 
wavelet frames. This can potentially be a problem because it might be difficult or even 
impossible to find a dual frame with wavelet structure. Since we exhibit pairs of dual 
wavelet frames, this issue is avoided. 

The principal importance of having a dual generator of the form (1) is that it will 
inherit properties from ip preserved by dilation and linearity, e.g. vanishing moments, 
good time localization and regularity properties. For a more complete account of such 
matters we refer to [15]. 

In the rest of this introduction we review basic definitions. A frame for a separable 
Hilbert space % is a countable collection of vectors {fj}jej for which there are constants 
< G\ < C 2 < oo such that 

Ci ll/ll 2 < 52\(f, fjf < C 2 H/ll 2 for all fen. 

jeu 

If the upper bound holds in the above inequality, then {fj} is said to be a Bessel 
sequence with Bessel constant C%- For a Bessel sequence {fj} we define the frame 
operator by 

s-.n^n, sj ^ (/../})./}. 

jeU 

This operator is bounded, invertible, and positive. A frame {fj} is said to be tight if 
we can choose C\ = C 2 , this is equivalent to S = C\I where I is the identity operator. 
Two Bessel sequences {fj} and {gj} are said to be dual frames if 

jeu 

It can be shown that two such Bessel sequences are indeed frames. Given a frame {fj}, 
at least one dual always exists; it is called the canonical dual and is given by {S~ 1 fj}. 
Only a frame, which is not a basis, has several duals. 

For / G L 1 (R n ) the Fourier transform is defined by /(£) = J Rn f(x)e~ 2m ^' x 'dx with 
the usual extension to L 2 (R n ). 

Sets in R™ are, in general, considered equal if they are equal up to sets of measure 
zero. The boundary of a set E is denoted by BE, the interior by E°, and the closure 
by E. Let B G GL n (M). A multiplicative tiling set E for {B^ : j 6 Z} is a subset of 
positive measure such that 

(2) R n \[j B j (E) = and \B j (E) n B\E) \ =0 for I ^ j. 

jez 

In this case we say that {B^(E) : j G Z} is an almost everywhere partition of R n , or 
that it tiles R™. A multiplicative tiling set E is bounded if E is a bounded set and ^ E. 
By i?-dilative periodicity of a function /: R n — > C we understand f(x) = f(Bx) for 
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a.e. x G R n , and by a I?-dilative partition of unity we understand Yljezf(^ x ) = 1; 
note that the functions in the "partition of unity" are not assumed to be non-negative, 
but can take any real or complex value. 

A (full-rank) lattice T in R n is a point set of the form T = PZ n for some P G GL n (R). 
The determinant of T is d(T) = |det P\; note that the generating matrix P is not unique, 
and that d(T) is independent of the particular choice of P. 



2. The general form of the construction procedure 

Fix the dimension n£N. We let A G GL n (M.) be expansive, i.e., all eigenvalues of 
A have absolute value greater than one, and denote the transpose matrix by B = A 1 . 
For any such dilation A, we want to construct a pair of functions that generate dual 
wavelet frames for some translation lattice. Our construction is based on the following 
result which is a consequence of the characterizing equations for dual wavelet frames 
by Chui, Czaja, Maggioni, and Weiss [6, Theorem 4]. 

Theorem 2.1. Let A G GL n (M) be expansive, let T be a lattice in R n ; and let \& = 
{■01, . . . , i^l}, ^ — {ipi, ■ ■ ■ i^Pl} C L 2 (R n ). Suppose that the two wavelet systems 
{D AJ T y ipi : j G Z, 7 G T, I = 1, . . . , L} and {D A iT 7 ijji : j G Z, y G V, I = 1, . . . , L} form 
Bessel families. Then {D^jTjijji} and {D^jT^ipi} will be dual frames if the following 
conditions hold 



( 3 ) ^L.'EhiB' = d(T) a.e. £ G R n , 

1=1 jez 

L A 

(4) + 7) = a. e. £ G R" for 7 G T* \ {0}. 
i=i 

Proof. By £ = _B J 'u; for j G Z, condition (4) becomes 

L ^ 

(5) J2^i( Bju) )^l( Biu} + 7) = a - e - w £ K" for 7 G T* \ {0}. 

We use the notation as in [6], thus A(A,T) = {a 6 I" : 3(j,7) G Z x T* : a = B~ j y} 
and JA,r(a) = G Z x T* : a = S"^}. Since J A r(«) C Z x (r* \ {0}) for any 

a G A(A, T) \ {0}, equation (5) yields 



^2MB j u)i)t(Bi(u + B-iy)) = a.e. u G R" 



d(T) 

for a^O. By 1^(0) = Z x {0}, we can rewrite (3) as 

l 



^$i{B j u)$i{Bi{u + B-iy)) = 1 a.e. u G R n , 



V ' (J,7)6/A,r(0) i=l 
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using that B~ 3 j = for all j G Z. Gathering the two equations displayed above yields 

i L . 

^7n S E^ 5Jw ^( 5J ( u + B "^) = ^ a.e. u; G W 1 , 

for all a G A (A, T). The conclusion follows now from [6, Theorem 4]. □ 

The following result, Lemma 2.2, gives a sufficient condition for a wavelet system to 
form a Bessel sequence; it is an extension of [3, Theorem 11.2.3] from L 2 (IR) to L 2 (W l ). 

Lemma 2.2. Let A G GL n (R) be expansive, T a lattice inW 1 , andcf) G L 2 (R"). Suppose 
that, for some set M cR n satisfying \J l&z B l (M) = R n , 

(6) C ^^i su PEE I t + 7) < oo. 

Then the wavelet system {Dj^T^cj)} 3 -gz 7e r a Bessel sequence with bound C 2 - Further, 
if also 

(7) C'i = 4^i^|E|^)| 2 -E E 1^(^00(^ + 7) | ] > 0, 

1 ' \j&L jez 7 er*\{o} / 

holds, then {Dj^T^} j&^&v ^ s a frame for L 2 (M. n ) with frame bounds C\ and C%. 

Proof. The statement follows directly by applying Theorem 3.1 in [5] on generalized shift 
invariant systems to wavelet systems. In the general result for generalized shift invariant 
systems [5, Theorem 3.1], the supremum/infimum is taken over M™, but because of the 
5-dilative periodicity of the series in (6) and (7) for wavelet systems, it suffices to take 
the supremum/infimum over a set Mcl" that has the property that U l&z B l (M) = R n 
up to sets of measure zero. □ 

Theorem 2.1 and Lemma 2.2 are all we need to prove the following result on pairs of 
dual wavelet frames. 

Theorem 2.3. Let A G GL n (M) be expansive and ip G L 2 (R n ). Suppose that ip is a 
bounded, real-valued function with suppip C U^ =0 B~ J (E) for some d G No and some 
bounded multiplicative tiling set E for {B^ : j G Z} ; and that 

(8) E *P( Bj = 1 f° r fl.e. £ G 

Let bj G C for j = —d, . . . , d and let m = max {j : bj ^ 0} and m = — min {j : bj ^ 0}. 
Take a lattice V in W n such that 

d m+d 

(9) (\jB~ j (E)+J) D |J B~ j (E) = <t) for all 7 G T* \ {0}, 

j=o 3=-m 
and define the function <p by 

m 

(10) (j){x) = d{T) ^ b j l det A \' j ^{A~ j x) for x G R n . 

j=-m 
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If bo = 1 and bj + b-j = 2 for j = 1,2, ... ,d, then the functions ip and generate dual 
frames {£> AJ T 7 ^}j e z, 7 er and {A^^LeZ/rer for L 2 (R n ). 

Proof. On the Fourier side, the definition in (10) becomes 

771 

4>(0 = d(T) £ b^{B^). 

j=-m 

Since ip by assumption is compactly supported in a "ringlike" structure bounded away 
from the origin, this will also be the case for 0. This property implies that ip and will 
generate wavelet Bessel sequences. The details are as follows. The support of ip and 



is 



supp^ C \jB- j (E), 

j=0 



m+d 



supp C |J B-' J (E). 

j=-m 

Note that < m, m < d. The sets {B^(E) : j 6 Z} tiles W 1 , whereby we see that 



(12) 
and, 
(13) 



suppVK^'-) nB' d (E) 



for j < and j > d, 



supp 4>(B j -) fl B d (E) = for j < —m and j > m + d. 



Since m, m > 0, condition (9) implies that ■\p{B^) , il){B^ + 7) = for j > and 
7 G T* \ {0}. Therefore, using (12), we find that 



Y, 1^(^)^(^ + 7)1 = (^0) 2 < 00 for £ G IT d (£). 

jez 7 er* i=o 

An application of Lemma 2.2 with M = B~ d (E) shows that ip generates a Bessel 
sequence. Similar calculations using (13) will show that generates a Bessel sequence; 
in this case the sum over 7 G T* will be finite, but it will in general have more than one 
nonzero term. 

To conclude that ip and generate dual wavelet frames we will show that condi- 
tions (3) and (4) in Theorem 2.1 hold. By 5-dilation periodicity of the sum in condition 
(3), it is sufficient to verify this condition on B~ d (E). For £ G B~ d (E) we have by (12), 



d(T) 



3=0 



HO [bom + bMBO + ■■■ + bdHB d 
+ 4(B£) \b_M0 + b i>(B0 + ■■■ + b d _MB d 



+ ■ 



+ ^(B d b-MO + ■■■ + b^MB^O + b ^(B d O 
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and further, by an expansion of these terms, 

d 



3,1=0 
d 



3=0 3,1=0 
3>l 



Using that bo = 1 and bj_\ + = 2 for j ^ I and j, I = 0, . . . , d, we arrive at 

c( d 



V 7 J6Z i=0 j,l=0 

3>l 



a— n / \„^C7 / 



i=o 7 v jez 

exhibiting that ip and satisfy condition (3). 

By (11) we see that condition (9) implies that the functions </> and + 7) will have 
disjoint support for 7 G V* \ {0}, hence (4) is satisfied. □ 

Remark 1 . The use of the parameters bj in the definition of the dual generator together 
with the condition b-j + bj = 2 was first seen in the work of Christensen and Kim [4] 
on pairs of dual Gabor frames. 

We can restate Theorem 2.3 for wavelet systems with standard translation lattice Z n 
and dilation A = P~ X AP, where P G GLJR) is so that T = PIT 1 . The result follows 
directly by an application of the relations D^Dp = DpD^j for j G Z and DpTpk = T^Dp 
for k G Z ra , and the fact that Dp is unitary as an operator on L 2 (IR n ). 

Corollary 2.4. Suppose ip, {bj}, A and T are as in Theorem 2.3. Let P G GL n (M.) be 
such that T = PU 1 , and let A = P~ l AP. Then the functions ip = Dpip and <p = Dp(f), 
where <j) is defined in (10), generate dual frames {DjjTkip} jez,fcez n an d {D^T^} j£% } kez n 
forL 2 (R n ). 

The following Example 1 is an application of Theorem 2.3 in L 2 (K 2 ) for the quincunx 
matrix. In particular, we construct a partition of unity of the form (8) for the quincunx 
matrix. 

Example 1. The quincunx matrix is defined as 



and its action on M 2 corresponds to a counter clockwise rotation of 45 degrees and a 
dilation by V^hx2- Define the tent shaped, piecewise linear function g by 
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FIGURE 1. Sketch of the triangular domains Jj, i - 



1,2,3,4,5. 



'-l + 2xi + 2x 2 , 


for (xi,x 2 ) 


e Ji, 


2x 2 , 


for (xi,x 2 ) 


e J 2 , 


< 2xi, 


for (xi,x 2 ) 


e J 3 , 


< 2-2xi, 


for (xi,x 2 ) 


e J 4 , 


2 - 2x 2 , 


for (xi,x 2 ) 


e J 5 , 




V 


otherwise, 





^(£i,k) 



where the sets J; are the triangular domains sketched in Figure 1. Note that the value 
at "the top of the tent" is g(l/2, 1/2) = 1. Define ip as a mirroring of g in the Xi axis 
and the X2 axis: 

f<?(6,6) for (6,6) e [0,oo) x [0,oo), 

5(6, -6) for (6, 6) e [0, oo) x (-oo, 0) , 
<?(-6, 6) for (6, 6) e (-oo, 0) x [0, oo) , 

U(-Ci,-6) for (6,6) e (-oo,0) x (-oo,0). 
Since the transpose B of the quincunx matrix also corresponds to a rotation of 45 
degrees (but clockwise) and a dilation by a/2/ 2X 2, we see that X^jez^C^'O = ^ 

We are now ready to apply Theorem 2.3 with £7 = [—1, l] 2 \5~ 1 ([— 1, l] 2 ) = [—1, 1] 2 \ 
Ii and d = 2; the set E is the union of the domians J4 and J5 and their mirrored versions. 
We choose 6_ 2 = 5_i = and bi = b 2 = 2d(r), hence m = and m = 2. Therefore, 

d m+d 

\jB-i(E), |J b~\e) c[-i,i} 2 , 

that shows that we can take T* = 2Z 2 or T = 1/2Z 2 , since ([-1, l] 2 + 7) fl [-1, l] 2 = 
whenever 0^76 2Z 2 . Defining the dual generator according to (16) yields 

(14) <f>{x) = (l/4)^(x) + (1/4)V(A-V) + (l/8)^(A- 2 x); 
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using that d(T) = 1/4, and we remark that is a piecewise linear function since this is 
the case for if). The conclusion from Theorem 2.3 is that ip and <fi generate dual frames 
{D A jT k/2 ip}j t kez and {D A jT k/2 (f)}j,k& for L 2 (R 2 ). 

The frame bounds can be found using Lemma 2.2 since the series (6) and (7) are 
finite sums on E; for {D A jT k /2ip} one finds C\ = 4/3 and Ci = 4. 

When the result on constructing pairs of dual wavelet frames is written in the gen- 
erality of Theorem 2.3, it is not always clear how to choose the set E and the lattice 
T. In Example 1 we showed how this can be done for the quincunx dilation matrix and 
constructed a pair of dual frame wavelets. In Section 3 and Theorem 3.3 we specify how 
to choose E and T for general dilations. The issue of exhibiting functions ip satisfying 
the condition (8) is addressed in Section 4. 

In one dimension, however, it is straightforward to make good choices of E and V as 
is seen by the following corollary of Theorem 2.3. The corollary unifies the construction 
procedures in Theorem 2 and Proposition 1 from [15] in a general procedure. 

Corollary 2.5. Let d G No, a > 1, and ip G L 2 (IR). Suppose that ip is a bounded, 
real-valued function with supp^ C [— a c , — a c-d_1 ] U [a c ~ d_1 , a c ] for some c G Z, and 
that 

(15) ^^( a ^) = l for a.e. £ G R. 

Let bj G C for j = —d, . . . , d, let m = — min {j : {bj ^ 0}}, and define the function <fi 
by 

d 

(16) 0(^) = bjO~^iJj(a~^x) for x G KL 

j=-m 

Let b G (0, a _c (l + a" 1 ) -1 ]. Ifbo — b and bj + b-j = 2b for j = 1,2, ... ,d, then ip and (j) 
generate dual frames {D a3 T bk ip} jMIj and {D^T bk (f)} jtk& for L 2 (R). 

Proof. In Theorem 2.3 for n = 1 and A = a we take E = [—a c , — a c_1 ] U [a c_1 , a c ] as the 
multiplicative tiling set for {a J : j G Z}. The assumption on the support of ip becomes 

d 

supp^ C [J a~ j {E) = \-a c , -a^ 1 ] U [a" - *" 1 , a c ]. 

3=0 

Moreover, since 

d 2d 

|J a-\E) C [-a c , a c ] , (J a~ J '(£) C [-a c+m , a c+m ] , 

jr'=0 j=-m 

and 

([-a c , a c ] + 7) n [-a c+m , a c+m ] = for | 7 | > a c + a c+m = a c (l + a m ), 

the choice V* = 6~ 1 Z for 6 _1 > a c (l + a m ) satisfies equation (9). This corresponds to 
r = bZ for < b < a~ c (l + a™)" 1 . □ 
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The assumptions in Corollary 2.5 imply that m G {0, 1, . . . , d}\ we note that in case 
m — 0, the corollary reduces to [15, Theorem 2]. 

3. A SPECIAL CASE OF THE CONSTRUCTION PROCEDURE 

We aim for a more automated construction procedure than what we have from The- 
orem 2.3, in particular, we therefore need to deal with good ways of choosing E and 
T. The basic idea in this automation process will be to choose £ as a dilation of the 
difference between and B^ 1 ^*), where I* is the unit ball in a norm in which the 
matrix B = A f is expanding "in all directions"; we will make this statement precise in 
Section 3.1. This idea is instrumental in the proof of Theorem 3.3. 

3.1. Some results on expansive matrices. We need the following well-known equiv- 
alent conditions for a (non-singular) matrix being expansive. 

Proposition 3.1. For B G C7L n (R) the following assertions are equivalent: 
(i) B is expansive, i.e., all eigenvalues Aj of B satisfy |Aj| > 1. 
(ii) For any norm \ - \ on M. n there are constants A > 1 and c > 1 such that 

\B j x\ > l/c\ j \x\ for all j G N , 

for any x G M n . 

(Hi) There is a Hermitian norm | • ^ on M n and a constant A > 1 such that 

\B j x\* > X j \x\^ for all j G N , 

for any x G M n . 

(iv) 8 C X8 C BE for some ellipsoid 8 = {x G M n : \Px\ < 1} ; P G GL n (R), and 
A > 1. 

By Proposition 3.1 we have that for a given expansive matrix B, there exists a scalar 
product with the induced norm | • ^ so that 

\Bxl>X\xl for s 6 1", 

I I * — I I * " 

holds for some A > 1. We say that | - 1^ is a norm associated with the expansive matrix 
B. Note that such a norm is not unique; we will follow the construction as in the 
proof of [2, Lemma 2.2], so let c and A be as in (ii) in Proposition 3.1 for the standard 
Euclidean norm with 1 < A < |Aj| for i = 1, . . . ,n, where A« are the eigenvalues of 
B. For k G N satisfying k > 2 In c/ In A we introduce the symmetric, positive definite 
matrix K G GL n (R): 

(17) K — I + (B^fB' 1 + ■ • • + (B- k ) l B- k . 

The scalar product associated with B is then defined by (x, y)^ = x l Ky. It might not 
be effortless to estimate c and A for some given B, but it is obvious that we just need to 
pick k G N such that B t KB — X 2 K becomes positive semi-definite for some A > 1 since 
this corresponds to (KBx, Bx) > A 2 (Kx,x), that is, \Bx\l > A 2 \x\l for all x G M n . 

We let I* denote the unit ball in the Hermitian norm | • I,, = |i^ 1//2 -| associated with 
B, i.e., 

(18) /, = {i 6 1" : \xl < 1} = {x G W 1 : \K 1/2 x\ < l} = {x G M n : x l Kx < 1} , 
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and we let O* denote the annulus 

O* = h \B-\h). 

The ringlike structure of is guaranteed by the fact that B is expanding in all direc- 
tions in the | • ^ norm, i.e., 

(19) J* C \h C B(Q, A > 1, 

which is (iv) in Proposition 3.1. We note that by an orthogonal substitution I* takes 
the form {i e I" : ^ix\ + • • • + < 1}, where /ij are the positive eigenvalues of K 
and x = Qx with the ith column of Q G 0(n) comprising of the ith eigenvector of 
K. The annulus O* is a bounded multiplicative tiling set for {B^ : j G Z}. This is a 
consequence of the following result. 

Lemma 3.2. Let B G GL n (M.) be an expansive matrix. For i ^ there is a unique 
j G Z so that B^x G O*; that is, 

(20) R n \ {0} = [J £ J '(0*) c%om£ tfraon. 

Proof. From equation (19) we know that {B l (I*)}i e z is a nested sequence of subsets of 
R n , thus 

B l (Q\B l - 1 (Q = B l (0*), /GZ, 

are disjoint sets. Since \B~^x\ Jf < A~ J \x\^ and l-B-'x^ > A J \x\^ for j > and A > 1, we 
also have 
i 

|J 5 m (a) = \ = {x G R n : \B~ l x\* < 1 and > 1} 

m=-l+l 

DjiGl": A - ' |a?|„ < 1 and A' « > 1} = {x G M" : A"' < « < A'} . 
Taking the limit / — > oo we get (20). □ 
Example 2. Let the following dilation matrix be given 

'3 -3 N 



(21) A 







Here we are interested in the transpose matrix B = A 1 with eigenvalues fi\^ = 3/2 ± 
i>/3/2, hence £ is an expansive matrix with l/xi^l = v3 > 1. The dilation matrix B is 
not expanding in the standard norm | • \ 2 in IR ra , i.e., I2 B(I 2 ), as shown by Figure 2. 
In order to have B expanding the unit ball we need to use the Hermitian norm from (iii) 
in Proposition 3.1 associated with B. In (17) we take k = 2 so that the real, symmetric, 
positive definite matrix K is 

K = I + (B-ifB-i + (B-y B -* = (28/9 W 9 ) , 

and let (x, y)^ := x t Ky. The choice k = 2 suffices since it makes B l KB — X 2 K semi- 
positive definite for A = 1.03 and thus 

\Bxi>X\xL, xeM 2 , 
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Figure 2. Boundaries of the sets I 2 , B(I 2 ), B 2 (I 2 ), and B 3 (I 2 ) marked 
by solid, long dashed, dashed, and dotted lines, respectively. Note that 
I 2 \ B(I 2 ) is non-empty, and even I 2 \ B 2 {I 2 ) is non-empty. 

holds for A = 1.03. 

Figure 3 and 4 illustrate that B indeed expands the Hermitian norm unit ball /* in 
all directions. We also remark that the Hermitian norm with k = 1 will not make the 




FIGURE 3. The unit ball in the Hermitian norm | • 1^ associated with 
B and its dilations -B(/*), B 2 (I*), B 3 (I*). Only the boundaries are 
marked. 



dilation matrix B expanding in IR ra ; in this case we have a situation similar to Figure 2. 
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Figure 4. A zoom of Figure 3. Boundaries of the sets I*, B(I*), B 2 (P), 
and -B 3 (/*) marked by solid, long dashed, dashed, and dotted lines, re- 
spectively. 



3.2. A crude lattice choice. Let us consider the setup in Theorem 2.3 with the set 
E = B c {0*) for some c G Z, where the norm | ■ 1^ = l-ft' 1 / 2 -] is associated with B. Let /i 
be the smallest eigenvalue of K such that i = y/l/Ji is the largest semi-principal axis 
of the ellipsoid i.e., i = ma.x x£ j t \x\ 2 . Then we can take any lattice T = PZ n , where 
P is a non-singular matrix satisfying 

(22) ||P|| 2 < 



£\\A%(1 + \\A^\\ 2 ) 



as our translation lattice in Theorem 2.3. To see this, recall that we are looking for a 
lattice T* such that, for 7 G T* \ {0}, 

(23) supp 4> H supp ?/>(• ± 7) = 0. 

For our choice of E we find that supp</> C £> c+ — (/*) and supp-?/' C B C (E). Since 
|P c+ -d 9 < ||5 C+ ^|L M 2 < ||S C+ -|U for any x G /*, 



and similar for B c x, we have the situation in (23) whenever |7| 2 > £(||A C || 2 + ||A C+ — | 
Here we have used that for the 2-norm ||A|| 2 = ||-B|| 2 . For z G Z™ we have 



2) 



\z\ 2 <\\p t \\ 2 \(P t )- i z\ 2 = \\p\\ 2 \(P t r i z\ 2 , 

therefore, by \z\ 2 > 1 for z 7^ 0, we have 

\( pt )~ lz \ 2 ^ nW~ for^GZ\{0}. 

Ik II2 
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Now, by assuming that P satisfies (22), we have 



bi 2 = {(py'zu > i/n > t (i + ii^y > i(\\a% + \\a c+ ^\\ 2 ) 

for ^ 7 = (P 1 )-^ G (P t )- 1 Z n = r*, hence the claim follows. 



A lattice choice based on (22) can be rather crude, and produces consequently a 
wavelet system with unnecessarily many translates. From equation (22) it is obvious 
that any lattice F = PZ n with ||P|| sufficiently small will work as translation lattice for 
our pair of generators ip and <fi. Hence, the challenging part is to find a sparse translation 
lattice whereby we understand a lattice V with large determinant d(T) := |detP|. In 
the dual lattice system this corresponds to a dense lattice V* with small volume d(T*) 
of the fundamental parallelotope /p* since d{T)d{T*) = 1. In Theorem 3.3 in the next 
section we make a better choice of the translation lattice compared to what we have 
from (22). 

Using a crude lattice approach as above, we can easily transform the translation 
lattice to the integer lattice if we allow multiple generators. We pick a matrix P that 
satisfies condition (22) and whose inverse is integer valued, i.e., Q := P^ 1 G GL n (Z). 
The conclusion from Theorem 2.3 is that {DAjTQ-i k ip} j^z,k& n an d {AtJ^Q- 1 ^} j<=z,kez n 
are dual frames. The order of the quotient group Q~ lr L n /IT 1 is |detQ|, so let {di : 
i = 1, . . . , |det Q\} denote a complete set of representatives of the quotient group, and 
define 



Since {D A jT Q -i k il)}j &jk& n = {D AJ T k ip} jGliMZ n^ eq , and likewise for the dual frame, the 
statement follows. 

3.3. A concrete version of Theorem 2.3. We list some standing assumptions and 
conventions for this section. 

General setup. We assume A G GL n (K) is expansive. Let | • | # = ( - , ■ )*' 2 be a Her- 
mitian norm as in (iii) in Proposition 3.1 associated with B = A f , let denote the unit 
ball in the | • l^-norm, and let K G GL n (M.) be the symmetric, positive definite matrix 
such that (x, y)^ = y t Kx. Let A := diag(Ai, . . . , A n ), where {Aj} are the eigenvalues of 
K, and let Q G 0(n) be such that the spectral decomposition of K is Q l KQ = A. 

The following result is a special case of Theorem 2.3, where we, in particular, specify 
how to choose the translation lattice T. Since we in Theorem 3.3 define T, it allows for 
a more automated construction procedure. 

Theorem 3.3. Let A, I*, K, Q, A be as in the general setup. Letd G M andip G L 2 (R n ). 
Suppose that %p is a bounded, real-valued function with supp^ C B C (IA \ B c ~ d ~ l (I*) for 
some ceZ, and that (8) holds. Take T = (1 /2)A c Q\/~K r L n . Then the function ijj and 
the function <p defined by 



L j=o J 

generate dual frames {D^jT^^j^z^&r and {Dj^T^^j^^^ for L 2 (M. n ) 
Remark 2. Note that d(T) = 2~ n |det A\ c (X l ■ ■ ■ A n ) 1/2 and \fX = diag( v / AT, . . . , v%)- 



* = {T d ^ : i = 1, . . . , |det Q\} , $ = {T di : i = 1, . . . , |det Q\} . 



(24) 




for x G R n , 
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Proof. The annulus O* is a bounded multiplicative tiling set for the dilations {B^ : j G Z} 
by Lemma 3.2, hence this is also the case for B c (0*) for c G Z. The support of ijj is 
supp^ C B C (I*) \ B c ~ d ~ l (I*) = \Jj_ B°~i(O*). Therefore we can apply Theorem 2.3 
with E = B c (0*), bj = 2 and b_j = for j = 1, . . . , d so that m = and m = d. The 
only thing left to justify is the choice of the translation lattice 17 We need to show that 
condition (9) with m = and fn = d in Theorem 2.3 is satisfied by T* = 2B c QA~ 1 l 2 TI l . 
By the orthogonal substitution x = Qx the quadratic form x l Kx of equation (18) 
reduces to 

Xixl + • • • + A n x^, 
where A« > 0, hence in the x, = Q t x coordinates I* is given by 

which is an ellipsoid with semi axes —3=, . . . , -4=. Therefore, in the x coordinates, 

(I* + 7) n /* = for 0^76 2A- 1 / 2 Z n , 
or, in the x coordinates, 

(/, + 7) n h = for ^ 7 G 2QA- 1 / 2 Z n . 
By applying B c to this relation it becomes 

(25) (B C (Q + 7) n £ C (J*) = for ^ 7 G V* = 2B c QA-^ 2 Z n , 

whereby we see that condition (9) is satisfied with m = and T* = 2B c QA~ 1/2 Z n . The 
dual lattice of r* is T = l/2A~ c QA 1 / 2 Z n . It follows from Theorem 2.3 that ip and <\> 
generate dual frames for this choice of the translation lattice. □ 

The frame bounds for the pair of dual frames {D^jT^if)} and {D^jT^^j^x.-yer 
in Theorem 3.3 can be given explicitly as 

Ci = -^r inf Xf^^) 2 ' C2 = 17^ sup E(^(^)) 2 , 
d(r) tEBc-*(o*) \ J diV) €eB c- d(0t) ^ V J 

and 

d = -L inf £ (#(B^0) 2 , C 2 = -L sup £ U(B j tj) 2 , 
tZ(r) teBc-*(o*)^ d \ ) d{T) &b°-*(o.) jt? d V 7 

respectively. The frame bounds do not depend on the specific structure of T, but only 
on the determinant of T; in particular, the condition number C2/C1 is independent of 

r. 

To verify these frame bounds, we note that equation (25) together with the fact 
supp^, supp0 C £> c (/*) imply that 

1>(Qj>(Z + 7) = kOkt + 7) = for a.e. (Gl" and 7 G T* \ {0}. 
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Therefore, by equations (12) and (13) with E = B c (0*), m = and m = d, we have 

E E 1^(^)^(^+7)1 = E l^'of = E (^o) 2 , 



7 er* jez j=o 



and 



E E |0(^'O0(^+7)| = E | = E (^^) 

iez 7 er* jez j=-d 

for £ G B c ~ d (0*). The stated frame bounds follow from Lemma 2.2. 

Example 3. Let A and K be as in Example 2. The eigenvalues of 7\~ are Ai = 
(26 + 2v / 65)/9 « 4.7 and A 2 = (26 - 2y/EE)/9 « 1.1. Let the normalized (in the 
standard norm) eigenvectors of K be columns of Q G 0(2) and A = diag(Ai, A2), hence 
Q l KQ = A. By the orthogonal transformation x = Qx the Hermitian norm unit ball 
7* becomes 

\ 2 / \ 2 

Xi \ fx 2 



which is an ellipse with semimajor axis 1/a/A~2 ~ 0.95 and semiminor axis l/y/Xi ~ 0.46. 
Since A" 1 / 2 = diag(l/\/Ai, I/VA2), we have 

(J* + 7) n /J = for ^ 7 G 2A- 1 / 2 Z 2 . 

By the orthogonal substitution back to x coordinates, we get 

|(/* + 7) n 7*| = for + 7 G 2QA- 1/2 Z 2 . 

Suppose that ijj is a bounded, real- valued function with supp?/> C B C (I*) \B c ~ d ~ 1 (I^) 
for c = 1 that satisfies the 5-dilative partition (8). Since c = 1 we need to take 
r* = 25 1 QA~ 1 / 2 Z 2 and Y = l/2A~ 1 QA 1 / 2 Z 2 , see Figure 5 and 6. 

3.4. An alternative lattice choice. Let the setup up and assumptions be as in The- 
orem 3.3, except for the lattice T which we want to choose differently. As in Section 3.2 
the dual lattice T* needs to satisfy (23) for 7 G T* \ {0}. We want to choose T* as dense 
as possible since this will make the translation lattice F as sparse as possible and the 
wavelet system with as few translates as possible. Since supp ip, supp <p C B C (I*), we are 
looking for lattices T* that packs the ellipsoids 7? c (7„,) + 7, 7 G F*, in a non-overlapping, 
optimal way. By the coordinate transformation x = A~ 1 ^ 2 Q t B~ c x, the ellipsoid B C (I*) 
turns into the standard unit ball I2 in M. n . This calculations are as follows. 

B C (Q = {B c x : \x\l < 1} = {x : \K 1,2 B- C x\\ < l} 



{x:\K l / 2 B- c B c QA-V 2 x\l<l} 

{x : <x,A- 1/2 Q*7i"QA- 1/2 £) 2 < l} = {x : \x\\ < l} 



and we arrive at a standard sphere packing problem with lattice arrangement of non- 
overlapping unit n-balls. The proportion of the Euclidean space R n filled by the balls is 
called the density of the arrangement, and it is this density we want as high as possible. 
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Figure 5. The dual lattice T* = 2B C QA~ 1 / 2 Z 2 for c = 1 is shown by 
dots, and the boundary of the set B C (I*) by a solid line. Boundaries 
of the set B C (I*) translated to several different 7 G F \ {0} are shown 
with dashed lines. Recall that supp i/j, supp 4> C -B c (i*), hence supp^fl 
supp^(- + 7) = for 7 6 T* \ {0}. 




Figure 6. The translation lattice T = (\j2)A c Qh v l 2 l? for c = 1. 



Taking T as in Theorem 3.3 corresponds to a square packing of the unit n-balls I 2 + k 
by the lattice 2Z n , i.e., k e 2Z n . The density of this packing is V n 2~ n , where V n is the 
volume of the n-ball: V 2n = n n /{n\) and V 2n +i = (2 2n+1 n\n n ) / (2n + 1)!. This is not 
the densest packing of balls in M n since there exists a lattice with density bigger than 
1.68n2 _n for each n / 1 [9]; a slight improvement of this lower bound was obtained 
in [1] for n > 5. Moreover, the densest lattice packing of hyperspheres is known up to 
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dimension 8, see [19]; it is precisely this dense lattice we want to use in place of 2Z n 
(at least whenever n < 8). 

In M 2 Lagrange proved that the hexagonal packing, where each ball touches 6 other 
balls in a hexagonal lattice, has the highest density n/y/T2. Hence using PI? with 

instead of 2Z 2 improves the packing by a factor of 



7T/V12 . r— . r 
^=4/Vl2 = 2/V3. 

It is easily seen that this factor equals the relation between the area of the fundamental 
parallelogram of the two lattices |det 2/2x2! / |det P\. In Figure 5 we see that each ellipse 
only touches 4 other ellipses corresponding to the square packing 2Z n ; in the optimal 
packing each ellipse touch 6 others. In R 3 Gauss proved that the highest density is 
7r / -\/T8 obtained by the hexagonal close and face-centered cubic packing; here each ball 
touches 12 other balls. 



4. Dilative partition of unity 

With Theorem 3.3 at hand the only issue left is to specify how to construct functions 
satisfying the partition of unity (8) for any given expansive matrix. In the two examples 
of this section we outline possible ways of achieving this. 

4.1. Constructing a partition of unity. As usual we fix the dimension n G N and the 
expansive matrix B G GL n (M.). In the examples in this section we construct functions 
satisfying the assumptions in Theorem 3.3, that is, a real- valued function g G L 2 (IR n ) 
with suppg C £> c (/*) \ B c_d_1 (l*) for some c G Z and d G N so that the 5-dilative 
partition 

(26) ^^"0 = 1 fora.e. £gM", 

holds. 

In the construction we will use that the radial coordinate of the surface of the ellipsoid 
(9-B J (/*), j G Z, can be parametrized by the n — 1 angular coordinates 81, ... , Q n -\- The 
radial coordinate expression will be of the form h(6i, . . . ,^ n -i) -1 ^ 2 f° r some positive, 
trigonometric function h, where h is bounded away from zero and infinity with the 
specific form of h depending on the dimension n and the length and orientation of the 
ellipsoid axes. 

We illustrate this with the following example in M 4 . We want to find the radial 
coordinate r of the ellipsoid 

{x G M 4 : [x x lhf + (x 2 /£ 2 ) 2 + (x 3 /£ 3 ) 2 + (x 4 /£ 4 ) 2 = 1} , U > 0, i = 1, 2, 3, 4, 

as a function the angular coordinates 9\, 82 and 83. We express x = (xi,x 2 , £3, £4) G M 4 
in the hyperspherical coordinates (^61,82,83) G {0} UR + x [0, tt] x [0, 7r] x [0, 2ir) as 
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follows: 



£i = rcos0i, x 2 = r sm0icos6 2 , 

£3 = r sin 9\ sin 2 cos 63, x^ = r sin 9\ sin 2 sin #3 . 

1, . . . , 4, in the expression above and factor out r 2 to obtain 



Then we substitute Xj, i 
r 2 f(8 1 ,6 2 ,6 3 ) = l, where 

(27) m,eM = £?ccx 2 e 1 

+ 

The conclusion is that r = r (81, 62,63 



% 2 sin 2 0! cos 2 # 2 



sin 2 6>i sin 2 8 2 cos 2 3 + £j 2 sin 2 8 X sin 2 2 sin 2 3 . 



f{@i, #2, 63 



-1/2 



Example 4. For d = 1 in Theorem 3.3 we want g G C^R") for any given s G NU{0}. 
The choice d — 1 will fix the "size" of the support of g so that suppg C B C (I*)\B C ~ 2 (I :¥ ) 
for some c G Z. Now let rj = ri(0i, . . . , n -i) and r 2 = r 2 (0i, . . . , 0„_i) denote the radial 
coordinates of the surface of the ellipsoids d.B c_1 (i*) and dB c (I*) parametrized by n — 1 
angular coordinates #1, . . . , n _i, respectively. 

Let / be a continuous function on the annulus S = B c (0*) satisfying f\aB c - 1 (i») = 1 
and f\dB c (h) — 0. Using the parametrizations r 1; r 2 of the surfaces of the two ellipsoids 
and fixing the n— 1 angular coordinates we realize that we only have to find a continuous 
function / : [r l5 r 2 ] — >■ R of one variable (the radial coordinate) satisfying /(ri) = 1 and 
/(r 2 ) = 0. For example the general function / G C°(S) of d variables can be any of the 
functions below: 

r 2 — r 



(28a) 


/(x) = /(r, 


01,.. 


• • , 0n-l) = 


(28b) 


f(x) = /(r, 


0i,. 


• • , 0n-l) = 


(28c) 


/(x) = /( r , 


0i,.. 


• • , 0n-l) = 


where r = 


|ar| G [ri,r 2 ], 


0i,. 


• • , 0n-2 G 


ri(0i,...,0 


n _i) and r 2 = r 2 (0i 


, • • • , 0n-l)- 


even belong 


58 to C\S). 






Define g 


G L 2 (R) by: 












f(Bx) 


(29) 






) 











r 2 - ri 
(r 2 — r) 2 
(r 2 - ri) 3 



(2(r-ri) +r 2 -n), 



r— ri 



h + o COS 7r( 
2 2 ^r2— ri- 



[0, 7r], and n _! G [0, 27r); recall that r\ = 
. In definitions (28b) and (28c) the function / 



for x G B C -\Q\B C - 2 (Q, 
iox xeB c (Q\B c -\h), 
otherwise. 

This way g becomes a 5-dilative partition of unity with suppg C B C (I^) \ B C ~ 2 (I*), so 
we can apply Theorem 3.3 with ip = g and d = 2. 

We can simplify the expressions for the radial coordinates r 1; r 2 of the surface of the 
ellipsoids 95 c_1 (/*) and dB c (I t .) from the previous example by a suitable coordinate 
change. The idea is to transform the ellipsoid 5 C_1 (/*) to the standard unit ball 7 2 by 



a first coordinate change x = A 1 ^ 2 Q t B 



-c+l 



x. 



This will transform the outer ellipsoid 



B C (I*) to another ellipsoid. A second and orthogonal coordinate transform x = Q)x 
will make the semiaxes of this new ellipsoid parallel to the coordinate axes, leaving 
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the standard unit ball J 2 unchanged. Here Q, comes from the spectral decomposition 
of A~ 1 B~ 1 , i.e., A~ 1 B~ 1 = Q'jA/Q,. In the x coordinates r% — 1 is a constant and 
r 2 = f^ 1 / 2 with / of the form (27) for n = 4 and likewise for n^4. 

In the construction in Example 4 we assumed that d = 1. The next example works 
for all d G N; moreover, the constructed function will belong to C^°(R n ). 

Example 5. For sufficiently small 5 > define Ai, A 2 C M. n by 

A 1 = B c - d - 1 (I*) + 5(0,5), 
A 2 + 5(0,5) = B c {h). 

This makes A2 \ Ai a subset of the annulus B C (I*) \ i? c_d_1 (/ !|t ); it is exactly the subset, 
where points less than 5 in distance from the boundary have been removed, or in other 
words 

A 2 \ Ax + 5(0, 5) = B C (Q \ B c - d -\Q. 

For this to hold, we of course need to take 5 > sufficiently small, e.g. such that 
Ax C rAx C A 2 holds for some r > 1. 

Let h G C£°(lR n ) satisfy supp h = 5(0,1), h > 0, and f hdfi = 1, and define 
h$ = 5~ d h(5~ 1 -). By convoluting the characteristic function on A 2 \ Ax with hg we 
obtain a smooth function living on the annulus 5 C (/*) \ 5 c_a! ~ 1 (/*). So let p G C7^°(K rt ) 
be defined by 

P = h s * Xa 2 \Ai, 

and note that suppp = 5 c (i*) \ 5 c ~ d ~ 1 (i*) since supp hs = B(0,S). Normalizing the 
function p in a proper way will give us the function g we are looking for. We will 
normalize p by the function w. 



w{x) = y^p(5 



For a fixed x G M n \ {0} this sum has either d or d+ 1 nonzero terms, and w is therefore 
bounded away from and 00: 

3c, C > : c < w(x) < C for all x G R n \ {0}, 

hence we can define a function g G C^°(lR n ) by 

(30) c/(x) = ^H- for s G R n \ {0}, and, c/(0) = 0. 

w(x) 

The function g will be an almost everywhere 5-dilative partition of unity as is seen by 
using the 5-dilative periodicity of w. 

E = E 414 = E ^ = -k Erf**) = i- 

^-^ w(B^x) wix) w(x) 

Since p is supported on the annulus 5 c (i*) \ B c ~ d ~ l (I*), we can simplify the definition 
in (30) to get rid of the infinite sum in the denominator; this gives us the following 
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expression 

d 

g(x) = p(x)/ P{B j x) for xeW n \ {0}. 
j=-d 

We can obtain a more explicit expression for p by the following approach. Let r\ = 
r\(6i, . . . , 9 n -i) and r 2 = ^2(^1, . . . , n -i) denote the radial coordinates of the surface 
of the ellipsoids 9.B c_d_1 (/*) and dB c (I*) parametrized by n — 1 angular coordinates 
61,..., n -i, respectively. Finally, let p e C£°(R n ) be defined by 

p(x) = T)(\x\ - ri)r](r 2 - \x\), with n = ri(0i, . . .,0 n -i) and r 2 = r 2 (0i, . . .,0 n -i) 

where 0i, . . . , n -\ can be found from x, and 




x > 0, 
x < 0. 
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